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I show how applying a symplectic Gram-Schmidt orthogonalization to the normalizer of a quantum 
code gives a different way of determining the code's logical operators. This approach may be more 
natural in the setting where we produce a quantum code from classical codes because the generator 
matrices of the classical codes form the normalizer of the resulting quantum code. This technique is 
particularly useful in determining the logical operators of an entanglement-assisted code produced 
from two classical binary codes or from one classical quaternary code. Finally, this approach gives 
additional formulas for computing the amount of entanglement that an entanglement-assisted code 
requires. 
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I. INTRODUCTION 

The logical operators of a quantum code are impor- 
tant for fault-tolerant quantum computation [l[ and in 
the simulation of the performance of stabilizer quantum 
error-correcting codes [2j. The typical method for de- 
termining the logical operators of a quantum code is to 
employ Gottesman's algorithm Q. This algorithm brings 
the stabilizer of a quantum code into a standard form and 
extracts the logical operators from this standard form. 

In this paper, I provide another technique for deter- 
mining the logical operators of a quantum code. This 
technique begins with the normalizer of the code and is 
perhaps more natural in the setting where we produce 
a CSS quantum code from two classical additive binary 
codes [3, Q or where we produce a CRSS quantum code 
from one classical additive code over GF(4) ((J. In both 
of these cases, we possess the full normalizer of the re- 
sulting quantum code because the generator matrices of 
the classical codes form its normalizer. Additionally, I 
show that this technique is useful in determining the log- 
ical operators of an entanglement-assisted code @, Q , a 
type of quantum code to which Gottesman's algorithm 
does not apply. The technique is again natural in the 
setting where we import classical codes because the gen- 
erator matrices form the full normalizer of the resulting 
entanglement-assisted code. This technique also gives 
another way of computing the amount of entanglement 
that a given entanglement-assisted code requires (in this 
sense, this paper is a sequel to the findings in Ref. [§]). 
Specifically, I give a set of formulas, different from those 
in Ref. [9], for computing the entanglement that a given 
entanglement-assisted code requires. The computation 
of the formulas in this paper may be more efficient than 
those in Ref. @ and the efficiency depends on the param- 



eters of a given code. The technique for computing the 
logical operators and the entanglement formulas in this 
paper should be a useful addition to the quantum code 
designer's toolbox. 



II. REVIEW OF THE SYMPLECTIC 
GRAM-SCHMIDT ORTHOGONALIZATION 
PROCEDURE 

I first review the symplectic Gram-Schmidt orthogo- 
nalization procedure (SGSOP) before proceeding to the 
main development of this paper. The SGSOP is equiv- 
alent to finding a standard symplectic basis for a set of 
vectors [Toj . It has made appearances in the quantum 
information literature in entanglement-assisted quantum 
coding 0, Q , in finding an entanglement measure for a 
stabilizer state [ill ], in the simulation of measurement 
outcomes of stabilizer codes |12j . and in the decomposi- 
tion of subsystem codes [l3| . 

We begin with a set of m Pauli generators g\ , . . . , g m 
that do not necessarily form a commuting set. Consider 
generator g%. There are two possible cases for this gen- 
erator: 

1. Generator g\ commutes with all other generators 
92, ■ ■ ■ , gm- Set it aside in a "set of processed gener- 
ators", relabel generators . . . , g m as respective 
generators g±, . . . , g m -i and repeat the SGSOP on 
these remaining generators. 

2. Generator g\ anticommutes with some other gen- 
erator gj where 2 < j < m. Relabel gj as gi and 
vice versa. Modify the other generators as follows: 
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where the function / is equal to zero if its two ar- 
guments commute and is equal to one if its two 
arguments anticommute. (Note that we are not 
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concerned with the overall phases of any of the gen- 
erators in this paper.) The new generators g%, . . . , 
g m then commute with g\ and gi. Place g\ and g2 
into the "set of processed generators," relabel gen- 
erators ^3, . . . , g m as respective generators g%, . . . , 
.9m-2, and repeat the SGSOP on these remaining 
generators. 

The result of the SGSOP is to produce a set of m pro- 
cessed generators where 2c of them form anticommuting 
pairs (each pair commuting with the generators in the 
other pairs) and / of them commute with themselves and 
the c pairs of generators. It is a reversible algorithm, 
meaning that we can recover the original set of gener- 
ators by performing all of the steps in reverse, and its 
complexity is O (nm 2 ) . This complexity is the same as 
the complexity of Gottesman's algorithm for bringing the 
generating set into a standard form and determining the 
logical operators from this standard form Q . 

III. GENERAL STABILIZER CODES 

We now consider applying the SGSOP to the normal- 
izer N (S) of a stabilizer code [3| with stabilizer S. The 
result is that the SGSOP gives a different way for deter- 
mining the logical operators of a given stabilizer code. 

Suppose that we have a set (S) of n — k Pauli genera- 
tors where each Pauli generator in (S) is an n-fold tensor 
product of Pauli matrices. Suppose furthermore that (S) 
generates an abelian group. Suppose we also have a gen- 
crating set (N (S)) of p elements where p > n — k and 
(N (S)) generates the normalizer N (S) of the stabilizer- 
s'. The relation S C N (S) holds for this case because 
the set of generators in (S) generates an abelian group. 

We can perform the SGSOP on the generating set 
(N(S)). The SGSOP divides (iV (S)> into two different 
sets: the logical generating set (Sl) and the stabilizer 
generating set (S) . We can write the latter generating set 
as (S) because it generates the same group that the gen- 
erating set from the previous paragraph generates. The 
logical generating set (Sl) consists of k anticommuting 
pairs of generators that correspond to the logical opera- 
tors for the stabilizer code S. The size of the generating 
set (N (S)) is then p = n - k + 2k = n + k. 

We can also think about the code in terms of binary 
vectors that represent it [1, [l4j]. Suppose we have the 
normalizer matrix: 

N=[N Z \ N x ], 

where N is the binary representation of the generators 
in (N(S)). We can consider the "symplectic product 
matrix" Qjy Q where 

tl N = N z Nl + N x Nl, 

and addition is binary. It is then straightforward to show 
by the method of proof of Theorem 1 in Ref. [§] that 

-rank (Ojv) = k. 



The idea of the proof is that the application of the SG- 
SOP to the normalizer matrix TV transforms the matrix 
f2jv to the following standard form: 

k n — k 

0J©0[o], (i) 

»=1 j=l 

where the small and large © correspond to the direct sum 
operation, J is the matrix 



and [0] is the one-element zero matrix. Each matrix J in 
the direct sum corresponds to a symplectic pair and has 
rank two. Each symplectic pair corresponds to exactly 
one logical operator in the code. Each matrix [0] has rank 
zero and corresponds to an ancilla qubit. See Ref. [9( for 
more details. 



A. CSS Codes 

One major class of quantum codes is the class of CSS 
codes [3, Q ■ These codes allow us to import two "dual- 
containing" classical codes for use as a quantum code. 
The resulting quantum code inherits the error-correcting 
properties and rate from the classical codes. 

Here, we show that the SGSOP provides a way to de- 
termine the logical operators for a CSS code. It con- 
structs the full code (stabilizer and logical operators) 
from the generator matrices of the two classical codes, 
rather than exploiting the parity check matrices. 

Suppose we have two classical codes C\ and C2 that 
satisfy Q C\. Suppose their respective generator 
matrices are G\ and G2, and their respective parity check 
matrices are Hi and H2- The conditions HiGj = and 
H2G2 = follow from the definition of a parity check 
matrix, and the condition HxH^ = follows from the 
condition C2 Q C\. 

The typical method for constructing a CSS quantum 
code is to build the following stabilizer matrix: 



" Hi 








H 2 



The code forms a valid stabilizer matrix because H\H^ — 
(this condition corresponds to the requirement for a sta- 
bilizer code to consist of a commuting set of generators) . 

The number of generators in the CSS quantum code 
is equal to 2n — k\ — k2, and it therefore has n — 
(2n — fci — k-i) = k\ + &2 — n logical qubits with 
2 (fci + &2 — n) logical operators. The size of the normal- 
izer is thus (2n — ki — A^) + 2 (ki + k% — n) = k\ + fo. 

There is another way to build the code by exploiting 
the generator matrices G\ and G2. The below matrix is 
the normalizer matrix of the code: 






Gi ' 


G 2 






(3) 
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The conditions HiGj = and H2GJ — guarantee that 
the above matrix is a valid normalizcr matrix. It repre- 
sents the full normalizer for the code because the rows 
are linearly independent and the number of generators in 
the normalizcr is k\ + k 2 . 

Suppose that we run the SGSOP on the rows of the 
normalizer matrix in ([3]). The procedure then produces 
k\ + fox — n anticommuting pairs and n — k\ + n — k2 
commuting generators. The k\ + k 2 — n anticommuting 
pairs correspond to 2 (k\ + k 2 — n) logical operators for 
the code, and the n — k\ + n — k% commuting genera- 
tors correspond to the stabilizer generators of the code. 
This technique for obtaining the logical operators and 
the stabilizer generators in the case of a CSS code is per- 
haps more natural than applying Gottesman's method in 
Ref. Q because we already know the matrices G\ and G2 
when importing classical codes. 

By the same method of proof as Corollary 1 in Ref. @ , 
the following formula also holds for any CSS code 



rank (GiGj) = h + k 2 - 



(4) 



because this quantity captures the amount of anticom- 
mutativity in the generators in (|3j). 



B. CRSS Quantum Codes 

A CRSS quantum code is one that we create by im- 
porting a classical additive code over GF(4) [6|. 

Suppose a GF(4) code has parity check matrix H and 
generator matrix G where HG T = 0. Additionally, sup- 
pose that the code is dual under the trace product over 
GF(4), i.e., the following condition holds: 

tv{HH f } = 0, 

where the dagger symbol f is the conjugate transpose of 
matrices over GF(4) and the null matrix on the right has 
dimension (n — k) x [n — k). Note that the above trace 
operation is different from the standard matrix trace, but 
is rather an elementwise computation of the trace prod- 
uct over GF(4) [|. 

We can then create a quantum check matrix: 



ZJH 



(5) 



where 7 is the isomorphism that takes an n-dimensional 
vector over GF(4) to a n-fold tensor product of Pauli 
operators @. 

It is straightforward to determine the normalizer of 
such a code because we can construct it from the gener- 
ator matrix G: 



luG* 
uG* 



The above matrix is a valid normalizer for the resulting 
quantum code because 



tr ■ 



f 


' luH ' 




' loG* ' 






uJH 




log* 





tr 



luH 



[cuG T luG 7 



= tr 



HG T oJHG T 



jHG 1 



HG T 








(7) 



(6) 



where the first equality is by definition of the conjugate 
transpose, the second equality holds because lulu = 1, 
lulu = u>, and lulu = Tu for elements of GF(4). The last 
equality follows from the condition HG T = 0. 

We then can perform the SGSOP on the above nor- 
malizer matrix in The result of the algorithm is 
to produce 2 (2k — n) anticommuting pairs and 2 (n — k) 
generators that commute with themselves and the pairs. 
This method again determines the logical operators for 
the imported code and may be a more natural alterna- 
tive to Gottesman's method because we already know the 
matrix G that determines the normalizer for the resulting 
CRSS quantum code. 

By the same method of proof as Corollary 2 in Ref. 0] , 
the following formula holds for a trace-orthogonal code 
imported from GF(4): 

rank (GG 1 ) =2k-n. 



IV. ENTANGLEMENT-ASSISTED QUANTUM 
CODES 



We can extend this procedure to an entanglement- 
assisted code as well. This alternative procedure for 
determining the logical operators is useful here because 
Gottesman's method does not apply to a nonabelian 
group of Pauli generators. 

Suppose that we have a set (S) of m Pauli generators 
that represents a generating set of operators with desir- 
able error-correcting properties. The generating set (S) 
does not necessarily generate an abelian group. Suppose 
we also have a generating set (N (S)) of p elements. This 
generating set (N (S)} generates the normalizer N (S) of 
the stabilizer S. In general, the relation S C N (S) does 
not necessarily hold for this case because the set of gen- 
erators in (S) is not an abelian group. 

We can perform the SGSOP on the generators in (S). 
This procedure divides {S) into two sets of generators: 
the entanglement generator set (Se) and the isotropic 
generator set (Si). The set (Se) contains anticommut- 
ing pairs of generators (where generators different pairs 
commute), and the set (Si) contains a commuting set of 
generators that furthermore commute with the genera- 
tors in (Se)- The size of the generating set (Se) is 2c, 
and the size of the set (Si) is i so that m — i + 2c. 

We can also perform the SGSOP on the generating 
set (N (S)). The SGSOP divides this generating set into 
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two sets: the logical generating set (Sl) and the isotropic 
generating set (Si). We can write this isotropic gener- 
ating set as (Si) because it generates the same group as 
the group generated by the isotropic generating set from 
the previous paragraph. The size of the generating set 
(N (S)) is then p = i + 21. The generators in the logical 
generating set (Sl) are the logical operators for the code. 

We can again consider this procedure in the binary vec- 
tor representation. Suppose the following matrix speci- 
fics the normalizer matrix: 

N=[N Z \ N x ]. 

Consider the matrix fijv where 

n N = N z N%+N x N%, 

and addition is binary. By the same method of proof as 
Theorem 1 in Ref. Q , we can show that 

-rank (Ojv) = I. 

A. CSS Entanglement-Assisted Codes 

We now consider extending the development in Sec- 
tion IIVI to the formulation of a CSS entanglement- 
assisted code @, [1] ■ This extension gives a useful way for 
determining the logical operators of a CSS entanglement- 
assisted quantum code. Additionally we determine a for- 
mula, different from that in Corollary 1 of Ref. [9j, for 
computing the amount of entanglement that a given CSS 
entanglement-assisted quantum code requires. 

Consider two arbitrary additive binary codes C\ and 
C2 with respective parity check matrices H\ and Hi and 
respective generator matrices G\ and Gi. We can form 
a quantum check matrix as follows: 



" Hi 








H 2 



(8) 



The above matrix is not a valid stabilizer matrix, but 
we can run the SGSOP on it to determine how to add 
entanglement in order to make this code a valid stabi- 
lizer code. After adding ebits to the code, it then be- 
comes an entanglement-assisted stabilizer code. A CSS 
entanglement-assisted code constructed in this way has 
In — k\ — ki generators (2c of which form anticommuting 
pairs and 2n — k\ — ki — 2c form the isotropic generating 
set), ki + k,2 — n + c logical qubits, and consumes c ebits 

0- 

The normalizer of the generators in ([5]) is again 








G 2 






and represents the full normalizer of the generators in ([5]) 
because the number of generators is k\ + ki and the rows 



are linearly independent. This time, when we run the al- 
gorithm, we get 2 (fci + ki — n + c) anticommuting pairs 
and n — k\ + n — ki — 2c commuting generators. The 
commuting generators form the basis for the isotropic 
generating set of the code and the anticommuting gener- 
ators correspond to the logical operators of the code. 

By the method of proof of Corollary 1 in Ref. 0] , we 
can show that the following relation holds: 

rank(GiG 2 P ) = k% +k 2 -n + c. (9) 

This formula gives us another way to determine the 
amount of ebits that a given entanglement-assisted code 
requires because the parameters n, ki, and ki are known 
parameters. 

The complexity of the computation in (0) is O (kxk 2 n). 
A similar formula from Ref. [9J computes the amount of 
entanglement that the code requires: 

rank (HiH%) = c. 

The complexity of the above formula is 
O ((11 — fci) (n — ki) n). Thus, the formula in Q 
may provide a speedup for low-rate codes. It is up to 
the quantum code designer to decide which formula to 
use depending on the parameters of the code. 

B. CRSS Entanglement- Assisted Codes 

The method for constructing a CRSS entanglement- 
assisted code is again to import an additive classical code 
over GF(4). Suppose the imported code has parity check 
matrix H and generator matrix G so that HG T = 0. 
The imported code does not have to be orthogonal with 
respect to the trace product. We then create a quantum 
check matrix: 




where 7 is the same isomorphism as in (|S|). 

The normalizer of the code is again the same as in ([5]) . 
The normalizer matrix is 




and it is a valid normalizer for the same reasons as in ((T)) . 

Performing the symplcctic Gram-Schmidt orthogonal- 
ization procedure on the normalizer matrix, we then 
get 2 (2k — n + c) anticommuting pairs and 2 (n — k — c) 
commuting generators. This method again determines 
the logical operators for the imported code. 

By the same method of proof as Corollary 2 in Ref. 0] , 
we can show that the following formula holds: 

rank (GG f ) = 2k -n + c. (10) 

The formula gives another way for determining the 
amount of entanglement that an entanglement-assisted 
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quantum code requires because the parameters k and n 
are known parameters. 

The complexity of the computation of rank (GG^) is 
O (fc 2 n) . A similar formula from Ref. 0] computes the 
amount of entanglement that the codes requires: 

rank (HH f ) = c. 

The complexity of the above formula is 0((n — k) 2 n). 
Thus, the formula in (jlOp may again provide a speedup 
for low-rate codes. 



V. CONCLUSION 

I have shown how the normalizer of a quantum code, 
whether stabilizer or entanglement-assisted, gives a use- 
ful way for determining its logical operators. It is nat- 
ural to assume full knowledge of the normalizer in the 
case where we import classical codes to produce quan- 
tum codes. Additionally, this development gives formu- 
las for computing the amount of entanglement that an 
entanglement-assisted code requires, and these formulas 
are an alternative to those given in Ref. Q . It is straight- 
forward to extend these methods and formulas to qudit 



codes and continuous-variable codes, by glancing at the 
results in Ref. @. 

Future work includes determining how to apply these 
findings in the setting of quantum convolutional codes 
[IH fl(| and entanglement-assisted quantum convolu- 
tional codes [It}. Todd Brun and I have made attempts 
at this problem by devising the algorithm in Section IV-D 
of Ref. 18] for the CSS case and the expansion technique 
and polynomial symplcctic Gram-Schmidt algorithm in 
respective Sections IV and V of Ref. [l|| for the CRSS 
case. It is more complicated in the convolutional setting 
because we would like to maintain the periodic structure 
of the quantum convolutional code while still maintain- 
ing the standard symplectic relations. An algorithm that 
achieves both tasks would be a boon to the theory of 
quantum convolutional coding. 
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